Hamiltonian chaos in a coupled BEC — optomechanical cavity system 
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We study a hybrid optomechanical system consisting of a Bose-Einstein condensate (BEC) trapped 
inside a single-mode optical cavity with a moving end-mirror. The intracavity light field has a dual 
role: it excites a momentum side-mode of the condensate, and acts as a nonlinear spring that couples 
the vibrating mirror to that collective density excitation. We present the dynamics in a regime where 
the intracavity optical field, the mirror, and the side-mode excitation all display bistable behavior. 
In this regime we find that the dynamics of the system exhibits Hamiltonian chaos for appropriate 
initial conditions. 

PACS numbers: 42.50Pq, 37.30-(-i, 85.85.-|-j, 42.65.-k 



I. INTRODUCTION 

Cavity optomechanics is currently the focus of exten- 
sive theoretical and experimental investigations and has 
witnessed spectacular advances in the last few years, see 
Ref. [l[ and references therein. A central paradigm of 
the field has been the cooling of one or several modes 
of vibration of a movable mirror in an optical cavity to 
their quantum mechanical ground state. Rapid progress 
has been made towards this goal using micro- and nano- 
fabricated mirrors [2:], membranes [j] and zippers 
Equally impressive has been the demonstration of op- 
tomechanics effects in situations where the oscillating 
mirror is replaced by other physical objects, most no- 
tably perhaps an ultracold atomic gas Q or a BEC 
These experiments are first steps in opening another fron- 
tier broached by optomechanics - the coherent coupling of 
mechanical elements with atomic systems 0; B ; oi' more 
generally the interfacing of solid state physics with AMO 
science. It is expected that these developments will lead 
to basic advances, for instance in the understanding of 
the quantum-classical interface, as well as to novel appli- 
cations related to measurement and sensing , coherent 
control 0, and quantum information processing [ll[ . 

The present paper considers theoretically a system at 
the boundary between cavity QED, nanoscience, and ul- 
tracold science - a BEC trapped inside a single-mode 
optical cavity with a moving mirror, see Fig. [TJ This 
system was previously discussed in the limit where the 
coupling between the BEC and the optical field is weak, 
the object of that study being the ensuing bistability of 
the BEC Mott insulator-superfluid transition [l2l|. In 
contrast, we now consider a situation where the BEC is 
strongly coupled to the intracavity field. 

It has been shown experimentally for the case of a cav- 
ity with fixed end-mirrors that the recoil resulting from 
the interaction between the ultracold atoms and the in- 
tracavity field results in the excitation of a condensate 
side-mode that is (in the simplest case) formally iden- 
tical with a mechanical oscillator Q. Extending these 
results to the case of a resonator with a moving end- 
mirror, and therefore coupling a "microscopic mirror" to 
a nanoscale device, is particularly interesting in view of 



FIG. 1: Bose-Einstein condensate trapped inside a Fabry- 
Perot cavity with a moving end-mirror. Due to photon recoil 
the intracavity light field induces a matter-wave momentum 
side-mode in the condensate that is formally analogous to 
a second, microscopic moving mirror coupled to the "true" 
moving mirror via the optical field. 



the high level of experimental control achievable in this 
system. As such it provides a test-bed for a number of 
investigations ranging from the use of a BEC as a quan- 
tum sensor to characterize and control the state of the 
moving mirror, j3] to (conversely) the use of the moving 
mirror to manipulate the condensate [isj , and to funda- 
mental studies of the entanglement between macroscopic 
and microscopic objects |14!] including the effects of de- 
coherence. 

At the simplest level the moving mirror and the con- 
densate side-mode can be approximated as two harmonic 
oscillators coupled by the intracavity optical field, which 
acts as a nonlinear spring. In general, all three compo- 
nents of the system can exhibit a multistable behavior, 
which normally reduces to a bistable behavior for small 
mirror displacements. In the bad cavity limit that we 
consider here, where the optical field decay rate is faster 
than the decay rates of the mirrors, this bistable behavior 
can result in rich mirror dynamics. In particular, in the 
limit where the damping of the mirrors can be completely 
neglected, the mirrors can undergo Hamiltonian chaos for 
appropriate initial conditions |15j . We find an energy- 
dependent order to chaos transition similar to that seen 
in the double pendulum and in the Henon-Heiles system 

The remainder of this paper is organized as follows: 
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Section II introduces our model, and discusses the anal- 
ogy between the condensate side-mode and a moving mir- 
ror driven by radiation pressure. We then eliminate adia- 
batically the intracavity field to derive coupled equations 
that describe the dynamics of the coupled condensate- 
moving mirror system in the absence of dissipation. Sec- 
tion III discusses the dynamics of the system, and shows 
the appearance of Hamiltonian chaos for appropriate ini- 
tial conditions. We show Poincare surfaces of section and 
spectra that illustrate that point. Finally, section IV is 
a summary and conclusion. 



II. THE MODEL 

We model the combined mirror-condensate-optical 
field system via the Hamiltonian 



H — Hr + Ha 



(1) 



where He describes the intracavity field and its coupling 
to the moving mirror. Ha describes the condensate and 
its coupling to the light field, and Hd accounts for the 
coupling of the various subsystems to thermal reservoirs 
and hence to dissipation. 

The Hamiltonian He is given explicitly by [13] 



He = —hAcCi^a — h^a^aq 



hij. 



m , »2 



ihri (a 



(2) 

where q is the dimensionless position operator for the 
quantized mode of vibration of the mirror of effective 
mass m under consideration and p its momentum op- 
erator, with commutation relation = i, a and 
are the bosonic annihilation and creation operators for 
the optical cavity field of frequency cUp, with [a, a^] = 1, 
^ ~ Ucy^h/mujm/ L is the optomechanical coupling rate, 
and Ac = ujp — luc is the cavity-pump detuning. Finally, 
the Hamiltonian —ihr]{d — a^) accounts for the external 
pumping of the cavity mode. 

We describe the motion of the BEC trapped in the op- 
tomechanical cavity with a simple one-dimensional model 
in which the atomic motion along the cavity axis is quan- 
tized. For large light-atom detunings we can adia- 
batically eliminate the internal excited state dynamics of 
the atoms and ignore the effects of spontaneous emission. 
Considering in addition low enough atomic densities that 
two-body interactions can be neglected we have 



2ma dx^ 



+ KUqci) acos^ kx ) ^ {x) dx 
(3) 



where 4' [x) is a bosonic field annihilation operator, ma 
is the atomic mass, Uq = g^/Aa is the far off-resonant 
vacuum Rabi frequency and k — cjp/c is the wave number 
of the light field. We consider in the following the strong 
coupling regime Ng^/ |Aq| ^ k, [18], where N is the aver- 
age number of atoms and k is the decay rate of the cavity, 
in which case the cavity field is strongly dependent on the 



collective density distribution of the BEC. More specif- 
ically the back-action of the BEC and the mirror can 
change the photon number inside the cavity significantly 
(but without changing its spatial structure). (Note that 
we ignore for now the effect of the harmonic trap con- 
fining the atoms inside the Fabry-Pcrot in Hamiltonian 
®-) 

The photon recoil associated with the absorption and 
stimulated emission of light by the BEC results in the 
generation of a symmetric momentum side-mode at ±2fifc 
p . To account for this effect we expand the field operator 
in Eq. ^ as 



if {x) ~ {cq + \/2cos{2kx)c-^ /Vi, 



(4) 



where cq and 62 are bosonic annihilation operators for 
atoms in the zero-momentum state and side-mode com- 
ponents, respectively. With this expansion the Hamilto- 
nian ([3]) simplifies to 



H' 



2 



(^c^co -I- c^C2 j H a'a (^c^C2 + c^^cq j 



+ 4:huJrclc2, 



(5) 



where the first and second terms account for the poten- 
tial energy and the optical coupling between the zero- 
momentum mode and the side-mode of the BEC respec- 
tively. The last term is the energy of the side-mode, of 
frequency iuJr = 2hk'^ /ma- 
in the absence of particle losses we have 



4 Co 



4^2 



(6) 



where the approximate equality holds if the zero- 
momentum component of the condensate is only weakly 
depleted by the coupling to the momentum side-mode. 
Treating then that component classically via cq and 



N, the coupling term in Eq. (O reduces to 



'2NhUQ.^ 
: a' a 



C2 



h^smd^dC 



where 



Ssm — j- 



msm4:UJr 



and 



L'^NU^UJr 



(7) 

(8) 
(9) 

(10) 



is the effective mass of the side-mode. 

This indicates that the side-mode is formally analogous 
to a mirror whose motion is driven by radiation pressure, 
in full analogy to the cavity end-mirror, with position 
operator Q and momentum operator P = i{c\ — C2)j\f2 
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satisfying the canonical commutation relation [Q, P] = 
i. We call this effective mirror the "sm- mirror" in the 
following. (We see below that the presence of dissipation 
weakens somewhat the analogy between the side-mode 
and a moving mirror.) 

A complete description of the system must include the 
effects of dissipation, on the optical field, the damping 
of the sm-mirror, and the mechanical damping of the 
end-mirror. These processes, which are accounted for by 
the Hamiltonian in Eq. ([T]), can be incorporated via 
standard quantum noise operators [l^. They result in 
the coupled quantum Langevin equations 



— = (lA + i^q — i^smQ — k) a + 7] + v^dj 



da 
H 
dq 
IE 
dp 
'dt 
dQ 
dt 
dp 

m 



-uj„iq + h^a)a - j„ip + Jb, 

AuJrP - JsmQ + flM, 

-AuJrQ - ft^smd^d - IsmP + /2M, 



(11) 



where A — Ac — NUq/2, din is the Markovian input noise 
of the cavity field, 7™ is the mechanical energy decay rate 
of the end-mirror and fs is the associated Brownian noise 
operator 20]. The harmonic trapping potential of the 
condensate, which we have ignored so far, couples the 
±2fi,/c-momentum side-modes to other modes, resulting 
in the damping of the sm-mirror Q . It is described by 
the terms involving jsun with fiM and /2M the asso- 
ciated noise operators, assumed to be Markovian. We 
remark that unlike for physical harmonic oscillators such 
as the end-mirror, the damping of the sm-mirror is not 
simply frictional: a damping term appears now both in 
the Langevin equation for the momentum P and for the 
position Q. 

Both the end-mirror as well as the sm-mirror are driven 
by the radiation pressure of the intracavity field, which in 
turn depends on the positions of these oscillators. Hence 
the microscopic and macroscopic mechanical oscillators 
are coupled through the optical field, which effectively 
acts as a nonlinear spring connecting them. 

We focus in the following on the classical dynamics of 
the two mirrors by treating their positions and momenta 
as classical variables. We also assume that the optical 
field decays at the fastest rate. This allows one to elimi- 
nate it adiabatically by setting its time derivative to zero, 
giving the steady-state intensity Ictsp, 



A 



4^,(l+72^/(4a;,)2) 



(12) 

where d — s- a is the classical limit. 

It is known that in the absence of a condensate this 
system can exhibit radiation pressure optical bistability 



[2lj. In the present case, with effectively two moving 
mirrors, one finds easily that the steady-state positions 
of both the end-mirror, {qs), and the sm-mirror, {Qs) 
have the same bistability properties as the optical field, 
via the relations 



Qs 



4'^r(l+7sm/(4^r 



(13) 



Examples of bistability curves are shown in Fig.[2ja) and 
(b). 



III. CLASSICAL DYNAMICS 

We assume as already indicated that the intracavity 
field follows adiabatically the motion of both the end- 
mirror and the sm-mirror, and that the retardation ef- 
fects resulting from optical damping are negligible [23 |. 
Treating in addition the motion of the mirrors classically 
and considering time scales short enough that their me- 
chanical damping can be ignored, their dynamics is gov- 
erned by the coupled equations of motion 



fq 
dt^ 



d^Q 
dt^ 



(4w,.) Q - 



(14) 
2.(15) 



It is easily verified that these equations can be derived 
from the effective classical Hamiltonian 



with 



T 
V 



H,s = T + V 



—p^ H -P^ 

9 + ^^<3 



(16) 



2 



arctan 



(a + - Ln.Q) /k] ■ (17) 



The Hamiltonian also describes the motion of 

a particle in the two-dimensional potential V{q,Q). 
V {q, Q) is generally a double- well potential, as illustrated 
in Fig. [2fc). The two local minima A and B correspond 
to the pairs of stable values {qA,QA\ and {<z_b,(3_b} of 
Fig. [2l^a,b), and the saddle point C corresponds to the 
unstable values {qc, Qc}- The existence of such a saddle 
point is the source of chaotic motion in many classical 
models, such as e.g. the Henon-Heiles model, with the 
appearance of Hamiltonian chaos determined by the en- 
ergy of the system [l^. We can similarly expect energy- 
dependent transitions between regular and chaotic mo- 
tion in the present case. 
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(a) 



(a) 



(b) 
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FIG. 2: (Color online. Bistable equilibrium position of (a) 
the end-mirror and (b) the sm-mirror as a function of the nor- 
malized dimensionless pump intensity 77^ /k^. The solid lines 
correspond to stable solutions and the dashed lines to unsta- 
ble ones. In this example the cavity length is i = 10^*m,the 
cavity-pump detuning Ac — 2Tr x 15MHz, and the mass and 
frequency of the movable mirror are 10-^g and 2tt x 19KHz, 
respectively. The wave length of the pump field is Ap — 780nm 
and the decay rate of the cavity is «: = 27r x 1.3MHz, and the 
decay of the sm-mirror is zero. The total atom number — 
1.2 X 10* and the vacuum Rabi frequency Uo = 2n x S.lKHz. 
(c) Two-dimensional contour plot of the potential V {q, Q) at 
"rf I >^ — 1.8, illustrating the existence of two local minima 
A (<3 = 1.7,0 = -3.2) and B {q ^ 9.6, Q = -18.1), with 
potential energies Ea = US.ShiUrn and Eb = ISTAhcorn, re- 



spectively, and of a saddle point C {q - 
Ec = 157 Mllo^. 



5.9,0 



-11.1), with 



The dynamics of the coupled-mirrors system takes 
place in the four-dimensional phase space {p,q,P,Q). 
In the absence of dissipation, energy conservation con- 
strains the system trajectories to a three-dimensional vol- 
ume of this four-dimensional space. The phase-space re- 
gion occupied by the system trajectories can therefore 



q 

(c) 



q 

(d) 



FIG. 3: (Color online) . Poincare sections (g, p) of the coupled- 
mirrors system for (a) an energy E = UOhuim close to the 
local minimum energy Eb, (b) an energy E — 157hLo m near 
but below the saddle-point energy Ec, (c) E = IGOfto;™ above 
Ec, and (d) E = ITOhuJm significantly above Ec- Same pa- 
rameters as in Fig. 2. 



be conveniently visualized by plotting the solutions of 
Eqs. (fT4ll and (fT5|l in a three-dimensional space, for in- 
stance (p, q, Q). When the motion of the mirrors is reg- 
ular, the trajectories are confined to closed tori in three- 
dimensional space that translate into regular closed or- 
bits on Poincare sections. When the motion becomes 
chaotic, these tori start to deform according to the 
Kolmogorov-Arnold-Moser theory [l^. For chaotic mo- 
tion the trajectories have no special geometric shape, and 
the Poincare sections turn into a "chaotic sea" without 
any distinctive structure. 

In the following we concentrate for concreteness on 
Poincare sections in the [p, q) plane for Q — Qb and 
P > 0. When solving Eqs. and numerically for 
random initial conditions corresponding to a total en- 
ergy slightly above Eb but much lower than the saddle- 
point energy Ec, we observe that the Poincare section 
is characterized by a series of regular elliptic orbits, see 
Fig. [3lja). These orbits can be divided into two families 
that surround two points (p > and p < 0) on the line 
q = qB- The oscillations of the end- mirror and the sm- 
mirror are in phase for p > and out of phase for p < 0. 
The observations imply that at low energies the motion 
of the coupled mirrors is quasi-periodic. We note that 
similar Poincare sections have previously been observed 
in the case of integrable two-dimensional potentials with 
anharmonicities up to terms of fourth order (23j , as well 
as for a double pendulum at low energy . In the effec- 
tive particle analogy, the particle is trapped in potential 
well B. 
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(a) 

q 

4r 
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(C) 

q 




FIG. 4: (a) g as a function of time for the initial condition 
p(0) = q{0) = P(0) = (5(0) = and the input intensity 
tj'^/k^ — 1.8. (b) The power spectrum of q{t), illustrating the 
quasi-periodic motion of the cavity end-mirror. The chaotic 
motion of the end-mirror is displayed by (c) and (d) which are 
for the same parameters as (a) and (b), except that 'rf' j = 2. 
All other parameters as in Fig. [2] 



The situation starts to change when the total energy 
of the system is close to E^- We then find that some of 
the out-of-phase elliptic orbits break into a series of ellip- 
tical islands, see Fig. [3{b). This is an indication that the 
closed torii start to be destroyed and the motion is tend- 
ing to chaos. When the energy of the system is slightly 
larger than the saddle-point energy, as shown in Fig.[3]Jc), 
all out-of-phase elliptic orbits dissolve into a chaotic sea, 
although some of the in-phase elliptic orbits still survive. 
Finally, when the energy is increased still further, those 
surviving orbits also degenerate into islands or dissolve 
in the chaotic sea, see Fig.[31[d). 

Figure S] shows the motion and the power spectrum of 
the end-mirror in both the quasi-periodic and the chaotic 
regimes. The series of Lorentzian-shaped peaks in the 
power spectrum of Fig. IH^b) is a clear signature of quasi- 
periodic motion, and should be contrasted to the broad, 
structureless spectrum of Fig. IH^d) , characteristic of the 
chaotic motion of qit) shown in Fig. |l{c). (The oscilla- 
tions of Q(t), not shown in Fig.[4l are of course likewise 
chaotic). In practice the chaotic dynamics of the sys- 
tem can be observed by detecting the light transmitted 
or reflected by the cavity. This light carries experimental 



(a) 




-8^ 

FIG. 5: The same parameters as those in Fig|4] but with 
decay terms, 7m = 7c = O.lii^m (a) The position of the mirror 
as a function of time, (b) The virtual positon of the sm-mirror 
as a function of time. 



signatures both of the end-mirror motion |25l . |26| and of 
the sm-mirror dynamics [6|. 

In currently experimentally realizable systems the me- 
chanical decay of the end-mirror is typically of the order 
of perhaps 7m ^ IHz. A soft harmonic trap for the con- 
densate, with trapping frequency 27r x IHz gives likewise 
a sm-mirror damping rate of the order of ~ IHz. 
In addition, one needs to account for the condensate 
lifetime, limited primarily by three-body collisions. For 
^^Rb atoms the three-body loss rate constant has been 
measured to be iiTs = 5.8 x lO-^Ocm^Hz [27*1. For an 
atomic density of n ~ 10 cm" -3 and iV - 1.2 X lO", the 
three-body loss rate is therefore Kj,rP'N ^ IHz. Hence 
the motion of the system is justifiably described in the 
framework of Hamiltonian chaos for time scales shorter 
than a second. 

For larger decay rates (or longer times), there is in prin- 
ciple the possibility that the system undergoes dissipative 
chaos. To determine whether this is the case we have car- 
ried out a linear stability analysis of Eqs. (|14p . properly 
generalized to include the effects of mirror damping. We 
have found numerically that the system is always stable 
in the bad cavity limit considered here. The situation 
is different in the case of a good (high finesse) cavity, 
where the intracavity field can no longer be adiabatically 
eliminated. In this case, the system can enter a regime 
of dissipative chaos, consistent with previous studies of 
chaos in optomechanical systems with a single moving 
mirror 25, 26.] . Figure O shows an example of mirror 
dynamics in the presence of larger mirror damping. It 
illustrates the disappearance of the chaotic motion and 
the approach to a stable steady state resulting from the 
damping guiding the system toward a a local energy min- 
imum. 
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IV. CONCLUSION AND OUTLOOK 

In conclusion, we have described the dynamics of a 
hybrid optomechanical system where a BEC is strongly 
coupled to the intracavity optical field of a Fabry-Perot 
cavity with a movable end-mirror. The photon recoil as- 
sociated with the virtual absorption and emission of light 
by the condensate results in the formation of an effective 
side-mode mirror that is driven just like the end-mirror 
by radiation pressure. In addition the light field couples 
nonlinearly this effective mirror and the end- mirror. We 
have shown that in the classical limit this optomechani- 
cal system displays radiation pressure induced bistability. 
In the bad cavity limit, and for times short enough that 
mirror damping can be safely ignored, the coupled dy- 
namics of the mirrors can be described by an effective 
double-well potential. For appropriate initial conditions 
the system is then characterized by the onset of Hamil- 



tonian chaos. 

Future work will extend this discussion to the good 
cavity limit, where the onset of dissipative chaos can also 
be expected. We will also consider a full quantum me- 
chanical description of the mirrors motion, and study 
the optically induced quantum correlations and entan- 
glement between the sm- mirror and the end- mirror. Ad- 
ditional goals include the quantum control of the mirror 
motion by the atoms, and conversely of the atoms by the 
optomechanical system. 
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